Practical approach to the sign problem at finite theta- vacuum angle 
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We investigate a way of circumventing the sign problem in lattice QCD simulations with a theta-vacuum 
term, using the PNJL model. We consider the reweighting method for the QCD Lagrangian after the J7a(1) 
transformation. In the Lagrangian, the P-odd mass term as a cause of the sign problem is minimized. In order 
to find out a good reference system in the reweighting method, we estimate the average reweighting factor by 
using the two-flavor PNJL model and eventually find a good reference system. 
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I. INTRODUCTION 

Phenomena based on strong interaction have shown that 
charge conjugation C, parity P and time reversal T are good 
symmetries of nature. This means that quantum chromody- 
namics (QCD) should respect any combinations of the dis- 
crete symmetries. Among the discrete symmetries, the CP 
symmetry is not necessarily respected in QCD due to the in- 
stanton solution [1,2]. The instanton solution allows the QCD 
Lagrangian Cqcd to have a (9-vacuum term. The resulting 
Lagrangian is described as 

Cqcd = q f (j v D u + m f )q f + ^F^F^ 

■i«i£« m F!,t (1) 
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in Euclidean spacetime, where F® is the field strength of 
gluon. The vacuum angle 9 is a periodic variable with period 
2ir. It was known to be an observable parameter [3]. The QCD 
Lagrangian is transformed as Cqcd(9) —> £qcd(—9) by the 
P transformation. Indicating that the P and CP symmetries 
are preserved only at 9 = and ±7r; note that 9 = —ir is iden- 
tical with 9 = 7t. In the 9 vacuum, therefore, we must consider 
the P and CP violating interaction parameterized by 9. The- 
oretically we can take any arbitrary value between —it and ir 
for 9. Nevertheless, it has been found from the measured neu- 
tron electric dipole moment [4] that \9\ < 10 -9 Q5|— tTD - Why is 
9 so small in zero temperature (T)? This long-standing puzzle 
is called the strong CP problem; see for example Ref. Q8[| for 
the review. 

Around the deconfinement transition at T = T4, there is 
a possibility that P-odd bubbles (metastable states) arise and 
thereby regions of nonzero 9 are generated [9]. Thus 9 can be- 
come a function depending on spacetime coordinates (t, x). If 
P-odd bubbles are really produced at T « Aqcd, P and CP 
symmetries can be violated locally in high-energy heavy-ion 
collisions or the early universe. This finite value of 9 could be 



a new source of large CP violation in the early universe and 
a crucial missing element for solving the puzzle of baryogen- 
esis. 

In the early stage of heavy-ion collision, the magnetic field 
is formed, and simultaneously the total number of particles 
plus antiparticles with right-handed helicity is deviated from 
that with left-handed helicity by the effective 9(t, x). In this 
situation, particles with right-handed helicity move opposite 
to antiparticles with right-handed helicity, and consequently 
an electromagnetic current is generated along the magnetic 
field. This is the so-called chiral magnetic effect IU0l - [l3ll . 
The chiral magnetic effect may explain the charge separa- 
tions observed in the recent STAR results |0. Hot QCD with 
nonzero 9 is thus quite interesting. 

For zero T and zero quark-number chemical potential (/i), 
some important properties are showed on P symmetry. Vafa 
and Witten proved for 9 = that the vacuum is unique and 
conserves P symmetry lfl5ll . This theorem does not preclude 
the existence of P-odd bubbles. At 6 = %, QCD possesses P 
symmetry as mentioned above, but it is spontaneously broken 
at low T Jl6l[T7il . The spontaneous violation of P symmetry 
is called the Dashen mechanism 11611 . Although the mecha- 
nism is a nonperturbative phenomenon, the first-principle lat- 
tice QCD (LQCD) is not applicable for finite 9 due to the 
sign problem. The mechanism at finite T and/or finite fi 
was then investigated with effective models such as the chiral 
perturbation theory 1H8I - I23I1 . the Nambu-Jona-Lasinio (NJL) 
model i24T - l27ll and the Polyakov-loop extended Nambu-Jona- 
Lasinio (PNJL) model j28M30h . LQCD has no sign problem 
at imaginary 9. Very recently the region was analyzed with 
LQCD and 9 dependence of the deconfinement transition tem- 
perature was investigated [31]. 

In the previous work 113011 ^ we proposed a way of minimiz- 
ing the sign problem on LQCD with finite 9. The proposal is 
as follows. For simplicity, we consider two-flavor QCD. The 
QCD Lagrangian (Q~|i is transformed into 
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with 



M (9) = lv D v + m cos (9/2) + mi~f 5 sin (9/2) (3) 
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by using the J7a(1) transformation 



q 



(4) 



where the quark field q — (q u ,qd) has been redefined by the 
new one q'. The determinant A4(0) satisfies 



detM{0) = [det M(-t 



(5) 



indicating that the sign problem is induced by the P-odd (9- 
odd) term, 771175 sin (0/2). The difficulty of the sign problem 
is minimized in (O, since the P-odd term with the light quark 
mass m is much smaller than the dynamical quark mass of 
order tIqcd- Actually, it was found that the P-even conden- 
sates a'f — (q'fq'f) is much larger than the P-odd condensates 
rj'f = (q'fij5q'f)- The P-even condensates little change even 
if the #-odd mass term is neglected. We then proposed the fol- 
lowing reweighting method. The vacuum expectation value of 
operator O is calculated by 



(6) 
(7) 

(8) 
(9) 



(O) = J VAOdetM(9)e- s <> 
= J VAO'dctM, ci {e)e- s o 



with the gluon part S g of the QCD action and 
O' 

R(9) 



R(6)0, 
detM(9) 



detM rc f(0) : 



where R(9) is the reweighting factor and det A4 Ie f(0) is the 
Fermion determinant of the reference theory that has no sign 
problem. The simplest candidate of the reference theory is the 
theory in which the 0-odd mass is neglected. We refer this 
reference theory to as reference A in this paper. As discussed 
in Ref. [30], reference A may be a good reference theory for 
small and intermediate 9, but not for large 9 near tt. In refer- 
ence A, the limit of 9 = n corresponds to the chiral limit that 
is hard for LQCD simulations to reach. 

The expectation value of R(9) in the reference theory is 
obtained by 



(R(9)) 



Z 



(10) 



where Z (Z rc [) is the partition function of the original (refer- 
ence) theory. The average reweighting factor (R(9)) is a good 
index for the reference theory to be good; the reference theory 
is good when (R(9)) = 1. 

In this paper, we estimate (R(9)) with the Polyakov-loop 
extended NJL (PNJL) model in order to find out a good ref- 
erence theory. Lagrangian (|2) does not depend on 9 in the 
chiral limit. For the realistic case closer to this limit rather 
than the pure gauge limit, an effect of finite 9 is characterized 
by the current quark mass. Hence, it is expected the finite 9 
effect is well described by effective models such as the PNJL 
model. As shown in Sec. [Till reference A is good only for 
small 9, so we propose the good reference theory that satisfies 



This paper is organized as follows. In Sec. [TT] we recapit- 
ulate the two-flavor PNJL model and show how to calculate 
the pion mass and (R(9)) for the case of finite 9. Numeri- 
cal results are shown in Sec. [HI] Section [IV] is devoted to a 
summary. 



II. MODEL SETTING 

The two-flavor PNJL Lagrangian with the ^-dependent 
anomaly term is obtained in Euclidean spacetime by 

3 

C = q(~f,yD v + m )q - Gi ^ [(qr a q) 2 + (qij5T a q) 2 ] 

a=0 

-8G 2 [e w detq R q L + e- w detq L q R ] + U(T,4>,$*), 

(11) 

where P„ = d v — ii v ^A\j \ a /2 with the Gell-Mann matrices 
A a . The current quark mass mo satisfies mo = m u = uid, and 
To and r a (a — 1, 2, 3) are the 2x2 unit and Pauli matrices in 
the flavor space, respectively. The parameter G\ denotes the 
coupling constant of the scalar and pseudoscalar-type four- 
quark interaction, while G2 stands for that of the Kobayashi- 
Maskawa-'t Hooft determinant interaction J2l[32[] where the 
matrix indices run in the flavor space. 

The gauge field is treated as a homo gen eous and static 
background field in the PNJL model JUS US The 
Polyakov-loop <P and its conjugate <P* are determined in the 
Euclidean space by 



<£=~tr c (Z), #* = itr c (L), 



(12) 



where L = exp(iA4/T) with A4/T = di&g(cf> r ,(j)g,<fib) in 
the Polyakov-gauge; note that A a is traceless and hence <\> r + 
4>g + 4>b — 0. Therefore we obtain 



(p 



5«* 

1 



J4>g 
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e i(<l>r+4>g)\ 

We use the Polyakov potential U of Ref. [37]: 



-i(f>g 



(13) 



u 



T 



a(T) 



+b(T) ln(l - 6<P<P* + 4(<? 3 + <T 3 ) - 3{<P<P*) 2 



with 



(T)=a„ + ai (f) + a 2 (fy, b(T)= 63(f) 



To" 



'T " 



(14) 



(15) 



The parameter set in U is fitted to LQCD data at finite T in the 
pure gauge limit. The parameters except To are summarized 
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in Table J] The Polyakov potential yields a first-order decon- 
finement phase transition at T = T in the pure gauge the- 
ory. The original value of To is 270 MeV determined from the 
pure gauge LQCD data, but the PNJL model with this value 
of To yields a larger value of the pseudocritical temperature 
T c of the deconfinement transition at zero chemical potential 
than T c « 173 ± 8 MeV predicted by full LQCD [40-42]. 
Therefore we rescale Tq to 212 MeV so as to reproduce the 
LQCD result. Under the Ua(1) transformation (|4), the quark- 



with 



do 


CJl 


ai 


63 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. [37], All parameters are dimensionless. 



antiquark condensates are transformed as 

a = qq = cos((9/2)cr' + sin(0/2)»/, (16) 

r\ = qij 5 q = -sin(0/2)<7' + cos(0/2)»/, (17) 

a, = qnq = cos(6»/2)a' i + sin(0/2)7r<, (18) 

tt 4 = qi l5Tl q = -sin(0/2X + cos(0/2)7r,;, (19) 

where the condensates {tr', rf \a\, n'j} are defined by the same 
form as {a, rj, a,i, 71^} but q is replaced by q' . The Lagrangian 
density is then rewritten with q' as 



where G± — G\ ± G 2 and 



m(9) = moCOs(0/2) + moi75sin(0/2). 



(22) 



Making the mean field approximation and the path integral 
over the quark field, one can obtain the thermodynamic po- 
tential fi (per volume) for finite T: 



n 



= u + u 



-Tin [l + 3&e- /3E± + 



3E± 

+Tln [1 + 3$*e-' 3E ± + 3$e~ 2pE± + e" 3 * 3 ^] 



(23) 



C ± 2Vf), 



E± = yp- 
C = M 2 + N 2 + A 2 
D = {MA + NP) 2 + (Ax Pf > 0, 
M = m o cos(0/2) - 2G+a', 
N = m o sin(0/2) - 2G-T)', 
A = -2G-a, P = -2G + w, 



A = VA-A, P = VP- P, 

U = G + (a' 2 +7r' 2 ) + G-( V ,2 +a' 2 ), 



(24) 
(25) 
(26) 
(27) 
(28) 
(29) 

(30) 
(31) 



where the momentum integral is regularized by the three- 
dimensional momentum cutoff A, Following Ref s. i25ll26ll . 
we introduce a parameter c as G\ = (1 — c)G+ and G2 = 
cG+, where < c < 0.5 and G+ > 0. The present model 
thus has four parameters of mo, A, G+ and c. Assuming 
mo = 5.5 MeV, we have determined A and G+ from the pion 
decay constant f n = 93 MeV and the pion mass = 138 
MeV at vacuum. Although c is an unknown parameter, we set 
c = 0.2 here, since it is known from the model analysis on the 
77 — 77' splitting that c w 0.2 is favorable [43]. 

For finite 8, parity is broken explicitly, so it is not a good 
quantum number anymore. Hence P-even and P-odd modes 
are mixed with each other for each meson. The "pion" mass 
M,,. is defined by the lowest pole mass of the inverse prop- 
agator in the isovector channel. It agrees with the ordinary 
pion mass when = 0. Under the random phase approxima- 
tion l44ll . the inverse propagator is described by 



det[l - 2GP(M^)] = 0, 



C = q'{ivD v + m(d))q' -Gt^2 [(9~'W) 2 + {q'il^q'f] 

a=Q 

-8G 2 [detq' R q' L + det^gJJ + U (20) 
= q'{luD v + m(9))q' - G+ [{q'q'f + (q'i l5 Tq>) 2 ] 

-G_ [{q'fq 1 ? + (?W) 2 ] +W, (2D 



where 



G 



n(u 2 ) = 




n ss {uj 2 ) n SP {uj 2 ) 
n PS {oj 2 ) n pp {uj 2 ) 



(32) 

(33) 
(34) 



with 



and 



n pp = 
n ss = 
n SP = 

h -- 
h{u 2 ) -- 

ftj> = 

= 



iN f Nj! - 2N f N c (q 2 - 4N 2 )I 2 (uj 2 ), (35) 
AN f N c h - 2N f N c {q 2 - AM 2 )I 2 {lo 2 ), (36) 



n 



PS = -8N f N c MNI 2 {uj 2 ), 
d 3 p l-f+(E p )-fe(E p ) 
(2tt) 3 2E p 

d 3 p l-f+(E p )-f#(E p ) 
A (2^) 3 



E p (uj 2 



(<2>* 



4E%) ' 



1 + 3(#* +<p e -P E p)e-P E p 
($ + 2$*e-P E »)e-P E p -+ 



-3/3E p 
. p—30E p 



(37) 
(38) 

(39) 



(40) 
(41) 
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In this form, we can set a! = %' = 0, since we do not consider 
the isospin chemical potential. 

Applying the saddle-point approximation to the path inte- 
gral in the partition function, one can get 



(R(d)) 



detii/ref _ 



/3V(0-0 Ie {) 



detH 



(42) 



where (3 = l/T, Q (f2 rc f ) is the thermodynamic potential of 
the original (reference) theory in the mean-field level, and H 
(H le f) is the Hessian matrix in the original (reference) theory 
defined by JUIH 



Hi 



d 2 n 
MM 



{4h} = {a', r]',ai,n'i, 



(43) 
(44) 



For later convenience, the average reweighting factor (R(0)) 
is divided into two factors Rh and Rq: 



{R{d)) = R H Rc 



with 



Rh = 



detif, 



rcf 



detH ' 
Ra = e-W- **). 



(45) 

(46) 
(47) 



For an JVj x N T lattice, the four-dimensional volume (3V is 
obtained by 



pv = 



1 



(48) 



Here we consider N X /N T — 4 as a typical example, following 
Refs. JUIH. 

We consider the following reference theory that has no sign 
problem: 

C = q'M + m Tei (9))q' - G+ [(q'q') 2 + (q'i^Tq') 2 ] 
-G- [{q'fq'f + (q'l l5 q' f] + U. (49) 

Here m re f (8) is #-even mass defined below. We consider three 
examples as m re f (9). 



A. The first example is reference A defined by 

m ie {(6) = mp,(6) = m cos(9/2). 



(50) 



In this case, the P-odd mass is simply neglected from 
the original Lagrangian (fJTJ. 

B. The second example is reference B defined by 

m le {(9) = m B (8) 

= m o cos(0/2) + - {m sin((9/2)} 2 . (51) 
a 

In this case, we have added the mQ-order correction 
due to the P-odd quark mass. Here a is a parameter 
with mass dimension, so we simply choose a = M„. 
The coefficient of the correction term is m^/M^ = 
0.129 Me V. 



(a) 
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\ <R> 
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\ 1 
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Original 




| Ref. A 



Fig. 1: 9 dependence of (a) the average reweighting factor and (b) 
M n at T = 100 MeV for the case of reference A. 



C. The third case is reference C defined by 

m rcl {9) = m c {9) 

= mo cos(9/2) + ^0-sm 2 (9/2), (52) 
v' 

This case also has the mg-order correction, but a is dif- 
ferent from reference B. The coefficient of the correc- 
tion term is m M 2 /M 2 , =0.114 MeV. 

Reference C is justified as follows. The pion mass M v (9) 
at finite 9 is estimated from the chiral Lagrangian as I123ll : 



M 2 A0) - ^Icos(0/2)| 

Jtt 



where cto is the chiral condensate at T = 9 = and the coef- 
ficient lj is evaluated by the 1/N C expansion as 



h 



JL 

2M 2 ' 



The right-hand side of (T53T > is reduced to 

m M 2 . 



M 2 (9) 



f 2 

■I 7T 



m |cos(f5/2)| 



M 2 , 



sm 2 (9/2) 



(54) 



(55) 



Equation (55[ supports (|52] >. 



m. NUMERICAL RESULTS 
A. Mean field approximation 

If some reference system satisfies the condition that 
(R(9)) ~ 1, one can say that the reference system is good. 
As a typical example of the condition, we consider the case 
of 0.5 < (R{9)) < 2. This condition seems to be the mini- 
mum requirement. The discussion made below is not changed 
qualitatively, even if one takes a stronger condition. 

First we consider reference A. Figure [TJa) shows 9 depen- 
dence of (R(9)) at T = 100 MeV. The solid line stands for 
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(a) 



<R> 




Fig. 2: 9 dependence of (a) the average reweighting factor and (b) 
M n at T = 100 MeV for the case of reference B. 



(a) 
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Fig. 3: 9 dependence of (a) the average reweighting factor and (b) 
Mtt at T = 100 MeV for the case of reference C. 



(R(9)), while the dashed (dotted) line corresponds to Rh 
(Rq). This temperature is lower than the chiral transition tem- 
perature in the original theory that is 212 MeV at 9 = and 
204 MeV at 6 = it. As 9 increases from zero, (R(9)} also 
increases and exceeds 2 at 9 « 1.2. Reference A is thus good 
for 9 < 1.2. The increase of (R(9)) stems from that of Rq 
that depends on T. This means that the reliable 9 region in 
which 0.5 < (R(9)) < 2 becomes large as T increases. 

Figure [Tib) shows 9 dependence of at T = 100 MeV. 
The solid (dashed) line denotes M n in the original (reference 
A) system. At 9 = 7r, M w is finite in the original system, but 
zero in reference A. As a consequence of this property, Rh 
and (R{9)) vanish at 9 = tt; see Fig. \V[a). This indicates that 
reference A breaks down at 9 = tt, independently of T. 

The same analysis is made for reference B in Fig. |2] M v 
in reference B well reproduces that in the original theory for 
any 9, and (R(9)) satisfies the condition 0.5 < (R{9)) < 2 
for all 9. Since Rh ~ 1 in the most region of 9, (R{9)) is 
governed by Rq. Around 9 = tt, Rh becomes small but 
still has a nonzero value because ^ even at 9 = tt in 
reference B. Therefore, the simple estimation for m re f (9) ( BTT l 
gives available reference. 

Finally we consider reference C through Fig. [3] M n in ref- 
erence C well simulates that in the original theory, and (R(9)) 
satisfies the condition 0.5 < (R{6)) < 2 for all 9. This result 



is better than that in reference B. Therefore we can think that 
reference C is a good reference system for any 9. This is true 
for any temperature larger than 100 MeV. 



B. Effect of mesonic fluctuation 

Beyond the mean field approximation, we estimate an ef- 
fect of dynamical pion fluctuations by modifying the thermo- 
dynamic potential to 



fi = fl 



MF 



Q 



DF, 



(56) 



where SImf is the thermodynamic potential d23l l with the 
mean-field level. ^7df is the potential due to dynamical pion 
fluctuations 14611 . 



(27T) 



(57) 



where E„ = y p 2 + M%, with M w determined by solving 

(El. 

Figure g] shows 9 dependence of (R{9)) at T = 100 MeV 
for the case of reference C. The solid and dashed lines corre- 
spond to results with and without dynamical pion fluctuations, 
respectively. The effect makes (R) a little smaller and hence 
the reference C becomes slightly worse. However, the modi- 
fication is small , indicating that (R) is well evaluated by the 
mean-field approximation 



IV. SUMMARY AND DISCUSSION 

We have investigated a way of circumventing the sign prob- 
lem in LQCD simulations with finite 9, using the PNJL model. 
We have considered the reweighting method for the trans- 
formed Lagrangian (HJ. In the Lagrangian, the sign problem 
is minimized, since the P-odd mass is much smaller than the 
dynamical quark mass of order Aqcb- Another key is which 
kind of reference system satisfies the condition (R(9)) s» 1. 




Fig. 4: 9 dependence of the average reweighting factor at T — 100 
MeV for the case of reference C. Solid and dashed lines correspond 
to the result with and without dynamical pion fluctuation respec- 
tively. 
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We have then estimated (R(6)) by using the two-flavor PNJL 
model and have found that reference C may be a good refer- 
ence system in the reweighting method. 

Since the present proposal is based on the model analysis, 
it is then not obvious whether the proposal really works in 
lattice simulations. Therefore the proposal should be directly 
tested by lattice simulations. A similar test was made for two- 
flavor QCD with finite quark chemical potential /u, ll46l l47ll 
where lattice simulations have the sign problem. The aver- 
age reweighting factor, i.e., the average phase factor was eval- 
uated by lattice simulations at fi/T < 1 for T around the 
critical temperature of the deconfinement transition |fi7ll . The 
PNJL model well reproduces the lattice result, when the dy- 



namical correction due to mesonic fluctuations is made to the 
mean-field model calculation 14611 . It is thus interesting that 
the present proposal is directly tested by lattice simulations. 
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